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Abstract We have investigated theoretically the effects of the charging en- 
ergy to the normal metal insulator-superconductor (NIS) tunnel junction 
used as a thermometer. We demonstrate by numerical calculations how the 
charging effects modify NIS thermometry and how the voltage-to-temperature 
response and the responsivity |dV/dT| of a current biased thermometer are 
affected. In addition, we show that the responsivity of the thermometer can 
be modulated with an additional gate electrode. The maximum responsivity 
is achieved when the Coulomb blockade is maximal, i.e. with a closed gate. 
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1 Introduction 

Tunnel junction thermometry with normal metal-insulator-superconductor 
(NIS) junctions was discussed a while backpQ, but has only recently been 
widely used in low-temperature thermal transport experiments [2,3,4,5. 
Advantages of NIS thermometry are e.g. low self-heating, good responsivity, 
ease of integration into the system under study, small size, existing high- 
frequency read-out schemes [5] and the fact that it can be used as local probe 
for temperature in nanostructures. Ideally, NIS junctions can be considered 
as primary thermometers, since their voltage-to-temperature (or current- 
to-temperature) calibration curves are determined by only two parameters: 
the superconducting gap A and the tunneling resistance Rt, which can be 
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determined quite accurately from the measured current-voltage (/-V) char- 
acteristics of the junction. However, charging effects[7] complicate this pic- 
ture for small area (capacitance) tunnel junctions connected to small islands 
(small self-capacitance). In this case charging effects influence the voltage(or 
current)-to-temperature response of the thermometers, if the charging en- 
ergy Ec — e 2 /(2C) is of the same order of magnitude or larger than tem- 
perature, i.e. Ec > T. This limit is easily obtainable in small devices, such 
as the recently demonstrated heat transistor [5] or the hybrid single electron 
transistor [9, lOj, where charging energies Ec ~ IK have been demonstrated. 
In this paper we investigate theoretically how Ec affects the current biased 
voltage-to-temperature responsivity of the SINIS thermometer and how it 
can be tuned with an additional gate located in the close proximity of the 
junctions. For simplicity, we only concentrate on the case, where single parti- 
cle tunneling is taken into account, and higher order two-particle (Andreev) 
processes [TT1[T2] are left out in the discussion. In addition to the charging ef- 
fects, we also discuss the effect of the non-ideal single particle current caused 
by the broadened density of states due to the finite life-time of quasiparticles 
in the superconductor [13], which causes deviations from the ideal behavior 
in thermometry. 

2 SINIS thermometry without charging effects 

The non-linear current-voltage (I-V) characteristics (Fig. [T] (a)) of a NIS 
tunnel junction can be used for thermometry (see e.g. |14j for a good re- 
view). In practice, thermometry is typically carried out by current biasing 
the junction at a constant current and measuring the voltage response of the 
thermometer, which is only a function of temperature. To obtain a larger 
signal, i.e. increased responsivity for thermometer, a structure containing 
two junctions in series (SINIS) can be used instead of a single NIS junction. 
Measured I-V curves for a typical SINIS thermometer with Ec << ksT 
at different bath temperatures are shown in Fig. Q] (a). It can be seen from 
Figs. Q] (a) and (b) that the voltage responsivity dV/dT of the thermometer 
can be adjusted to be optimal for a certain temperature range by adjust- 
ing the bias current: With low bias-currents (~ 10 pA) there is a higher 
sensitivity at low-temperatures, while with higher bias-currents (~ 100 pA) 
sensitivity is gained at high temperatures, but lost at low temperatures. The 
best results with thermometry are obtained by repeating experiments with 
a few different bias points for different temperature ranges. However, higher 
bias currents may cause significant heating at the lowest operating temper- 
atures due to power dissipation to the normal metal. The practical upper 
limit for the bias current I depends on the resistivity of the normal metal 
material, but is at most / corresponding to voltage Vsinis ~ 2zi/e, after 
which significant heating is induced by the junction itself[14 . Overheating 
is a critical issue at low temperatures, especially with structures containing 
small normal metal islands or otherwise very well thermally isolated samples, 
such as suspended nanowires [T51l 1 6j . In addition to heating, the bias current 
can provide self-cooling of the normal metal island [T7lll8j. Usually, SINIS 
junctions designed for thermometry have small junction areas and therefore 
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small self-cooling effects, that can be neglected in the analysis (Fig. Q] (a)). 
However, even for larger junctions, the self-cooling can be avoided by choos- 
ing a proper bias point that never yields voltages near the optimal cooling 
point Vsinis ~ 2A/e. 

Notice also from Fig. [1] (b) that typically the measured SINIS response 
as a function of the refrigerator (bath) temperature deviates from the theory 
calculation (Eq. © below) at T < 150 mK. This deviation is most likely 
due to noise heating of the electron gas, i.e. coupling of the electron gas into 
its electromagnetic environment causing overheating of the electrons so that 
-^electron 7^ Ibath- This noise heating power, typically ~ 10 fW is dependent 
on the filtering of the lines in the cryostat and the electrical impedance of 
the junctions. 
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Fig. 1 (Color online) (a) Typical measured sub-gap current -voltage characteristics 
of a SINIS probe junction at different temperatures (from right to left, lowest 
temperature on the right). Dashed horizontal lines from top to bottom correspond 
to bias currents of 100 pA and 10 pA, respectively, (b) SINIS thermometer voltage 
vs. bath temperature, black (bottom curve) and red (top curve) lines corrrespond 
to two different current bias points lOpA and 100 pA, respectively. Open circles 
represent the calculation from Eq. (J5J). 



The theoretical single particle current of a NIS junction derived from the 
tunneling Hamiltonian is given bv|19j 

I(V,T) = -^-J n s (E){f N (E-eV,T)-f s (E,T)]dE (1) 

where Rt is the tunneling resistance of the sample, ns{E) the density of 
states (DOS) of the superconductor and /jv,S the Fermi-Dirac distributions 
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in the normal metal and the superconductor, respectively. Interestingly, Eq. 
([1]) can be rewritten in a symmetric form 

/CT = rVf n s (E)[f N (E-eV,T)-f N (E + eV,T)]dE (2) 

where the Fermi-distribution of the superconductor is eliminated. Hence, 
the I—V is dependent only on the electron temperature of the normal metal, 
i.e. the NIS-junction probes directly the electron temperature of the normal 
metal and no additional knowledge of the temperature of the superconductor 
is needed. However, there is still an implicit dependence of the temperature 
of the superconductor, since the DOS depends on the superconducting gap, 
whose temperature dependence at T > 0.5Tc must be taken into account. 
Note also that if we relax the assumption of quasicquilibrium (Fermi-function 
distributions) for either terminal, this simplification is not valid anymore. 
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Fig. 2 (Color online) (a) I—V characteristics calculated at a temperature 50 mK 
with different values of the F parameter (smallest to largest value of FjA corre- 
sponds to curves from bottom to top). Values are shown in Fig. (b). Horizontal, 
dashed black line corresponds to the current bias that has been used for calculating 
the calibration curves in (b). This correponds to 10 pA for junction with Rt = 20 
kfi and A = 220 ueV. (b) SINIS thermometer voltage as a function of tempera- 
ture, corresponding to the same values of F as in (a) {r/A increases from top to 
bottom). 



Typically, real junctions show a non-exponential finite sub-gap current 
(cf. Fig. [1] (a) I ~ 1 pA) at low enough temperatures and voltages. This 
sub-gap current can also be modelled with Eqs. ((T])-([2|) by incorporating a 
broadened DOS, i.e. by taking into account the finite life-time of quasipar- 
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tides in the superconductor [13] ■ The broadened DOS is then written as 



where r describes the magnitude of broadening of the DOS. Typically for 
evaporated Al films we have F j A ~ 2 • 10~ 4 , consistent with values reported 
by other groups[20, 21j. FjA is known to depend strongly on the material 
quality and can be two orders of magnitude larger for evaporated Nb films 2 2J. 

Figures [2Ja) and (b) show a calculation based on Eq. @, how the broad- 
ening of the DOS affects the I-V characteristics and voltage-to-temperature 
response of a SINIS. It is clear that the broadening of the DOS starts to in- 
fluence the SINIS response mostly at low temperatures when FjA ~ 5 ■ 10 -4 , 
depending, of course, on the bias current value. At high biases, the broaden- 
ing can be neglected to rather high values of F but at lower biases the effect 
is stronger. The broadened DOS shows up in the Vsinis y S- T measurement 
qualitatively in the same way as noise heating, i.e. by leading to a saturation 
of the voltage at low temperatures. However, the two effects can be typi- 
cally resolved in the experiment because the value for r can be determined 
from the I-V characteristics, i.e. f changes the current deep in the sub-gap, 
whereas T changes the slope of the current rise at the gap edge (cf. Fig.[2ja)). 
Typically this broadening is not a problem for thermometry in evaporated 
Al films, because the bias point can be taken above the sub-gap current. 
However, if the SINIS thermometers are fabricated e.g. from evaporated Nb 
films where FjA may be relatively high, broadening can be a real problem 
reducing the responsivity of the thermometers. 

3 Tunneling current with charging energy 

Charging effects (Coulomb blockade) due to small capacitances of the tunnel 
junction have an effect on thermometry, especially in the limit Eq >> fcsT. 
We consider the case of two identical junctions in series (with the same 
capacitances C and tunneling resistances Rt), i- e. a symmetric SINIS ther- 
mometer. In addition, we discuss the case where a gate electrode is in close 
proximity with the normal metal island, so that the geometry is basically a 
hybrid sing le electron transistor (SET) [9] . The I-V characteristics of SINIS 
structures are naturally modified when the charging effects are taken into 
account, and can be controlled by applying a voltage V g to the gate [7]. The 
tunneling rates through a single junction with charging energy and biased 
with a voltage Vi can be written asfTU] 
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Here the E)^ = ±2E c {n + n g ± 0.5) ± eVi is the change in the electrostatic 
energy when an electron tunnels on to the island (+) and off the island (-) 
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through junction i = {L,R}, where L,R stand for left and right junction, 
respectively. Ec = e 2 /(2Cs)is the charging energy, where Cs ~ 2C + C g is 
the total capacitance of the island with C the junction capacitance and C g 
the gate capacitance, en is the excess quantized charge (n integer) on the 
island and en g = Q g the offset charge, which can be be varied continuously 
by the gate electrode voltage. 

The current through the device can be calculated using these tunneling 
rates by solving a Master equation with the detailed balance condition 



where T ± = r L - ± +T R ^ ± and P(n) is the occupation probability of the corre- 
sponding charge state n obeying the normalization condition X)^L-oo P( n ) = 
1. Once this equation is solved, the current through the island can be calcu- 
lated from the expression 



and I = II = Ir. It can be easily shown that Eq. reduces to Eq. j2]) in 
the limit Eq = 0. An interesting observation about the SINIS current with a 
charging energy included is that it cannot be written in a form that is fully 
independent of the superconductor temperature Ts, not even in the weak 
Coulomb blockade limit Ec < ksT. Hence, the conclusion of Eq. ([2|) where 
a NIS junction can be used as a normal metal electron thermometer without 
additional knowledge of Ts is only a special case when Ec « ksT. This 
effect may be important when considering high charging energy devices used 
for cooling [8] , where temperature differences between the normal metal and 
the superconductor occur because of the heat flow from the island to the 
leads. However, for typical temperature differences AT ~ 200 mK achieved 
for Al based SINIS coolers, the effect on the I—V characteristics is only visible 
in the sub-gap region, and with typical bias currents it can be neglected. 

4 Charging effects in thermometry 

In this section we present the main theoretical results of the charging effects 
on thermometry and discuss how the gate can modulate the responsivity of 
the thermometer. All results are presented in scaled units and in all calcula- 
tions the DOS broadening parameter T/A = 2- 10~ 4 has been used. Sub-gap 
I—V characteristics of a SINIS at temperature ksT/A = 0.02 with three dif- 
ferent values of charging energy Ec I A =0, 0.1 and 0.2 are shown in Fig. [3] 
It is clear that the charging energy effectively shifts the I-V curves to higher 
voltages, leading to a behavior that resembles an effective increase of the 
superconducting gap from A to maximally A + Ec (if n g = 0). Furthermore, 
since the DOS is non-zero within the gap, a weak fingerprint of the Coulomb 
staircase can be seen in the I—V curves in the sub-gap region at eV/A ~ 1.5. 

In order to study how thermometry is affected, we should investigate 
how the voltage-to-temperature response is modified. This is shown in figure 



r + (n)P(n) = r-(n+ l)P(n+ 1) 
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Fig. 3 (Color online) Calculated I-V characteristics of a SINIS thermometer with 
three different values of Ec/A — (black, solid), 0.1 (red, dashed) and 0.2 (blue, 
dot) and n g — 0. Two dashed horizontal lines correspond to the current bias values 
used in later calculations eR T I/A = 2.3 ■ 10" 3 and eR T I/A = 2.3 • 10" 2 . 



|4ja), calculated from Eqs. (J4j) — dTj) the same values of charging energy as 
used in calculating I-V characteristics for Fig. [3] The solid and dashed lines 
correspond to two different current bias values shown in Fig. [31 where the 
low bias eR T I/A = 2.3 ■ 1CT 3 , and the high bias eR T I/A = 2.3 • 10" 2 . These 
scaled current values correspond to 10 pA and 100 pA currents for junctions 
with Rt = 50 kfi and A — 220 peV, typical for thin film Al. With this value of 
A, the corresponding charging energies in Figs [3] and [4] are Ec/ks =0, 0.25 K 
and 0.5 K. The responsivities |dy/dT| of the SINIS thermometer calculated 
from the voltage vs. temperature curves are presented in Fig.[4jb). 

It can be seen from Fig. BJb) that the charging energy increases the 
responsivity of the thermometer, but it also changes the shape of the voltage 
vs. temperature curve: there appears a clear bump in the responsivity curve 
at around fceT = 0.15Z\ for the low bias (0.2Zi for the high bias) when the 
charging energy becomes appreciable. The shift of the effective gap is also 
clearly seen in Fig. 2] (a) in the zero temperature limit, which moves from 
eV — 2 A to eV = 2(A + Ec)- The responsivity curves also clearly show how 
the two different bias points have different optimal temperature ranges: at 
lower bias the thermometer has more responsivity for temperatures below 
the peak caused by charging effects (ksT/A < 0.15, which corresponds to 
T = 380 mK for the values of A and Rt used above), but quickly drops 
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Fig. 4 (Color online) (a) Calculated voltage response vs. temperature curves of a 
SINIS device for three different charging energies Ec/A = (black), O.f (red) and 
0.2 (blue), with n g = (maximal Coulomb blockade). Solid lines are with low bias 
{eR T I/A = 2.3 ■ 10~ 3 ) and dashed with high bias current (eR T I/A = 2.3 ■ f(T 2 ). 
(b) Responsivity |dV/dT| calculated from (a). 



to unpractically small values above that. However, at the higher bias the 
responsivity stays large up to much higher temperatures. 

Gate modulation of the voltage response vs. temperature curves and cor- 
responding responsivities are presented in Fig.[5Ja) and (b). Here the highest 
value of Eel A = 0.2 from the plot in Fig. 2] is used, and all curves are calcu- 
lated with the low current bias value eRrl/A = 2.3 • 10 -3 . Open and closed 
gate situations correspond to the top (magenta) and bottom curves (black) , 
respectively, and gate voltages between these two extrema are plotted with 
an interval of n g = 0.1. The open gate (n g — 0.5), interestingly, has the same 
zero temperature limit eV ~ 2A as a thermometer with zero charging energy, 
however the full temperature dependent responsivity is different. This can be 
seen most easily from the responsivity curves in Fig.[5jb), where we see that 
the responsivity bump caused by Ec still exist. Hence, the shape of the V 
vs. T curve is still different. Even more surprisingly, at low temperatures 
UbT < 0.02Z\ the responsivities with intermediate gate voltage values seem 
to merge, while the open and closed gate values are still well separated. This 
effect is more clearly shown in Fig|51 

Figure [5] shows the responsivity of the thermometer as a function of the 
gate charge number n g at different temperatures from ksT/A= 0.01 to 0.18. 
These curves were calculated with the same parameters as in Fig. [5J The 
low temperature results show sharply peaked values at the gate open and 
closed positions, while the intermediate values show flat regions, where the 
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Fig. 5 (Color online) (a) Calculated SINIS voltage vs. temperature curves for 
Ec/A = 0.2 and current bias eRrl/A = 2.3 • 10 -3 with different gate charge 
values varying from an open gate n 9 = 0.5(lowest, black curve)to a closed gate 
n g = (magenta, top curve) with interval of n g = 0.1. (b) Responsivities |dV/dT| 
calculated from (a) corresponding to the same gate values (lowest black curve open 
gate and top magenta curve closed gate). 



responsivity does not depend on the gate charge, as noted before. At higher 
temperatures, these sharp peaks broaden into a sinusoidal dependence whose 
amplitude decreases with increasing temperature. The largest change in the 
responsivity is obtained at the lowest temperatures, where responsivity can 
be enchanced by 30 % by closing the gate at ksT = 0.01Z\, for example. 



5 Relevance to experiments 

To decide in practice when one should consider Coulomb charging in SINIS 
thermometry, one needs to measure the value of the charging energy Eq and 
compare it with the temperature range of interest. The Eq measurement can 
easily be performed at 4.2 K in the weak Coulomb blockade (Eq < &bT) 
limit, where the size of the zero bias Coulomb blockade dip in the tunneling 
conductance spectrum AG depends on the charging energy by the relation 

m 

^ = (8) 
G T 3fc s T' v ; 

where Gt is the tunneling conductance around V — without the dip. 
As an example, figure [7] shows two measured Coulomb blockade dips at 4.2 
K with (a) a typical larger junction area (0.35 pm 2 ) Al/Cu/Al SINIS de- 
vice designed to act as a cooler [T5lfT6] . and (b) a smaller typical solitary 
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Al/Cu/Al SINIS thermometer with a junction area ~ 0.05 um 2 . From the 
measurements we obtain Ec/kB — 20 mK for the cooler sample, showing 
that it is in the limit of weak Coulomb blockade ksT > Ec for our exper- 
imentally achievable (dilution refrigerator) temperature range. The changes 
to the voltage-to-temperature response are within the experimental error in 
that case, and analysis can be carried out with the simpler theory of Eq. |(3J), 
without charging effects. 

However, the effects of charging energy are more observable in the second, 
smaller SINIS thermometer (5 um x 300 nm x 30 nm normal metal island), 
with a measured charging energy Ec/kB ~ 200 mK (Fig. [7] (b)). This type 
of thermometer is typically used e.g. for probing the bath or (local)phonon 
temperature during an experiment. The measured Ec corresponds to a value 
Eel A ~ 0.08 for Al, which leads to observable changes in the temperature 
response based on the calculations in the previous section. 



6 Conclusions 

We have studied the charging effects on SINIS tunnel junction thermometry 
and shown that for small enough junctions, the responsivity of the thermome- 
ter can be modulated with an additional gate electrode in close proximity to 
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Voltage (mV) Voltage (mV) 

Fig. 7 (Color online) The measured (black line) differential conductance spectrum 
of a large junction area SINIS device (0.35 pm 2 ) (a) and a smaller device (junction 

area 0.05 pi 2 ) (b) both at 4.2 K. Conductance is normalized with Gt = 1/Rt 
(red, dashed line). Solid (green/gray) line corresponds to a w eak-Coulomb blockade 
regime fit to the data in the lowest order in Ec /(foT) [14] . 

the junctions, with maximum responsivity achieved with the closed gate sit- 
uation. In addition, the shape of the voltage response vs. temperature curve 
changes when the charging energy is taken into account, leading to the con- 
clusion that in typical solitary SINIS thermometers charging effects must be 
taken into account in the conversion from measured voltage to temperature 
(calibration), and have to be computed numerically. The effect of an unknown 
offset charge will not lead to ambiguities in the analysis, as the shape of the 
voltage vs. temperature curve changes as a function of the offset charge, so 
that the value of the offset charge can be determined self-consistently. 
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